Bioinformatic methods in applied genomic research

Here I will focus on three main topics that best address and include the projects I have been working in during my three year PhD period that I have spent in different research laboratories addressing both computationally and practically important problems all related to modern molecular genomics. The first topic is the use of livestock species (pigs) as a model of obesity, a complex human dysfunction. My efforts here concern the detection and annotation of Single Nucleotide Polymorphisms. I developed a pipeline for mining human and porcine sequences. Starting from a set of human genes related with obesity the platform returns a list of annotated porcine SNPs extracted from a new set of potential obesity-genes. 565 of these SNPs were analyzed on an Illumina chip to test the involvement in obesity on a population composed by more than 500 pigs. Results will be discussed. All the computational analysis and experiments were done in collaboration with the Biocomputing group and Dr.Luca Fontanesi, respectively, under the direction of prof. Rita Casadio at the Bologna University, Italy. The second topic concerns developing a methodology, based on Factor Analysis, to simultaneously mine information from different levels of biological organization. With specific test cases we develop models of the complexity of the mRNA-miRNA molecular interaction in brain tumors measured indirectly by microarray and quantitative PCR. This work was done under the supervision of Prof. Christine Nardini, at the “CAS-MPG Partner Institute for Computational Biology” of Shangai, China (co-founded by the Max Planck Society and the Chinese Academy of Sciences jointly) The third topic concerns the development of a new method to overcome the variety of PCR technologies routinely adopted to characterize unknown flanking DNA regions of a viral integration locus of the human genome after clinical gene therapy. This new method is entirely based on next generation sequencing and it reduces the time required to detect insertion sites, decreasing the complexity of the procedure. This work was done in collaboration with the group of Dr. Manfred Schmidt at the Nationales Centrum für Tumorerkrankungen (Heidelberg, Germany) supervised by Dr. Annette Deichmann and Dr. Ali Nowrouzi. Furthermore I add as an Appendix the description of a R package for gene network reconstruction that I helped to develop for scientific usage (http://www.bioconductor.org/help/bioc-views/release/bioc/html/BUS.html).

Depicting the role of CDKN2A/ARF alterations in adult BCR-ABL1-positive acute lymphoblastic leukemia patients: from genomic deletions to prognostic impact

This 9p21 locus, encode for important proteins involved in cell cycle regulation and apoptosis containing the p16/CDKN2A (cyclin-dependent kinase inhibitor 2a) tumor suppressor gene and two other related genes, p14/ARF and p15/CDKN2B. This locus, is a major target of inactivation in the pathogenesis of a number of human tumors, both solid and haematologic, and is a frequent site of loss or deletion also in acute lymphoblastic leukemia (ALL) ranging from 18% to 45% 1. In order to explore, at high resolution, the frequency and size of alterations affecting this locus in adult BCR-ABL1-positive ALL and to investigate their prognostic value, 112 patients (101 de novo and 11 relapse cases) were analyzed by genome-wide single nucleotide polymorphisms arrays and gene candidate deep exon sequencing. Paired diagnosis-relapse samples were further available and analyzed for 19 (19%) cases. CDKN2A/ARF and CDKN2B genomic alterations were identified in 29% and 25% of newly diagnosed patients, respectively. Deletions were monoallelic in 72% of cases and in 43% the minimal overlapping region of the lost area spanned only the CDKN2A/2B gene locus. The analysis at the time of relapse showed an almost significant increase in the detection rate of CDKN2A/ARF loss (47%) compared to diagnosis (p = 0.06). Point mutations within the 9p21 locus were found at very low level with only a non-synonymous substition in the exon 2 of CDKN2A. Finally, correlation with clinical outcome showed that deletions of CDKN2A/B are significantly associated with poor outcome in terms of overall survival (p = 0.0206), disease free-survival (p = 0.0010) and cumulative incidence of relapse (p = 0.0014). The inactivation of 9p21 locus by genomic deletions is a frequent event in BCR-ABL1-positive ALL. Deletions are frequently acquired at the leukemia progression and work as a poor prognostic marker.

Analysis of the pig genome for the identification of genomic regions affecting production traits

The aim of this work was to identify markers associated with production traits in the pig genome using different approaches. We focused the attention on Italian Large White pig breed using Genome Wide Association Studies (GWAS) and applying a selective genotyping approach to increase the power of the analyses. Furthermore, we searched the pig genome using Next Generation Sequencing (NSG) Ion Torrent Technology to combine selective genotyping approach and deep sequencing for SNP discovery. Other two studies were carried on with a different approach. Allele frequency changes for SNPs affecting candidate genes and at Genome Wide level were analysed to identify selection signatures driven by selection program during the last 20 years. This approach confirmed that a great number of markers may affect production traits and that they are captured by the classical selection programs. GWAS revealed 123 significant or suggestively significant SNP associated with Back Fat Thickenss and 229 associated with Average Daily Gain. 16 Copy Number Variant Regions resulted more frequent in lean or fat pigs and showed that different copies of those region could have a limited impact on fat. These often appear to be involved in food intake and behavior, beside affecting genes involved in metabolic pathways and their expression. By combining NGS sequencing with selective genotyping approach, new variants where discovered and at least 54 are worth to be analysed in association studies. The study of groups of pigs undergone to stringent selection showed that allele frequency of some loci can drastically change if they are close to traits that are interesting for selection schemes. These approaches could be, in future, integrated in genomic selection plans.

Caratterizzazione genomica di virus influenzali suini in Italia mediante next generation sequencing

I sottotipi H1N1, H1N2 e H3N2 di influenza A virus sono largamente diffusi nella popolazione suina di tutto il mondo. Nel presente lavoro è stato sviluppato un protocollo di sequenziamento di c.d. nuova generazione, su piattaforma Ion Torrent PGM, idoneo per l’analisi di tutti i virus influenzali suini (SIV). Per valutare l’evoluzione molecolare dei SIV italiani, sono stati sequenziati ed analizzati mediante analisi genomica e filogenetica un totale di sessantadue ceppi di SIV appartenenti ai sottotipi H1N1, H1N2 e H3N2, isolati in Italia dal 1998 al 2014. Sono stati evidenziati in sei campioni due fenomeni di riassortimento: tutti i SIV H1N2 esaminati presentavano una neuraminidasi di derivazione umana, diversa da quella dei SIV H1N2 circolanti in Europa, inoltre l’emoagglutinina (HA) di due isolati H1N2 era originata dal riassortimento con un SIV H1N1 avian-like. L’analisi molecolare dell’HA ha permesso di rivelare un’inserzione di due amminoacidi in quattro SIV H1N1 pandemici e una delezione di due aminoacidi in quattro SIV H1N2, entrambe a livello del sito di legame con il recettore cellulare. E’ stata inoltre evidenziata un’elevata omologia di un SIV H1N1 con ceppi europei isolati negli anni ’80, suggerendo la possibile origine vaccinale di questo virus. E’ stato possibile, in aggiunta, applicare il nuovo protocollo sviluppato per sequenziare un virus influenzale aviare altamente patogeno trasmesso all’uomo, direttamente da campione biologico. La diversità genetica nei SIV esaminati in questo studio conferma l’importanza di un continuo monitoraggio della costellazione genomica dei virus influenzali nella popolazione suina.

Genomic and non genomic effects of elevated concentration of anabolic steroids in human neuronal cells

Nandrolone and other anabolic androgenic steroids (AAS) at elevated concentration can alter the expression and function of neurotransmitter systems and contribute to neuronal cell death. This effect can explain the behavioural changes, drug dependence and neuro degeneration observed in steroid abuser. Nandrolone treatment (10-8M–10-5M) caused a time- and concentration-dependent downregulation of mu opioid receptor (MOPr) transcripts in SH-SY5Y human neuroblastoma cells. This effect was prevented by the androgen receptor (AR) antagonist hydroxyflutamide. Receptor binding assays confirmed a decrease in MOPr of approximately 40% in nandrolonetreated cells. Treatment with actinomycin D (10-5M), a transcription inhibitor, revealed that nandrolone may regulate MOPr mRNA stability. In SH-SY5Y cells transfected with a human MOPr luciferase promoter/reporter construct, nandrolone did not alter the rate of gene transcription. These results suggest that nandrolone may regulate MOPr expression through post-transcriptional mechanisms requiring the AR. Cito-toxicity assays demonstrated a time- and concentration dependent decrease of cells viability in SH-SY5Y cells exposed to steroids (10-6M–10-4M). This toxic effects is independent of activation of AR and sigma-2 receptor. An increased of caspase-3 activity was observed in cells treated with Nandrolone 10-6M for 48h. Collectively, these data support the existence of two cellular mechanisms that might explain the neurological syndromes observed in steroids abuser.

Computational methods for genome screening

Motivation An actual issue of great interest, both under a theoretical and an applicative perspective, is the analysis of biological sequences for disclosing the information that they encode. The development of new technologies for genome sequencing in the last years, opened new fundamental problems since huge amounts of biological data still deserve an interpretation. Indeed, the sequencing is only the first step of the genome annotation process that consists in the assignment of biological information to each sequence. Hence given the large amount of available data, in silico methods became useful and necessary in order to extract relevant information from sequences. The availability of data from Genome Projects gave rise to new strategies for tackling the basic problems of computational biology such as the determination of the tridimensional structures of proteins, their biological function and their reciprocal interactions. Results The aim of this work has been the implementation of predictive methods that allow the extraction of information on the properties of genomes and proteins starting from the nucleotide and aminoacidic sequences, by taking advantage of the information provided by the comparison of the genome sequences from different species. In the first part of the work a comprehensive large scale genome comparison of 599 organisms is described. 2,6 million of sequences coming from 551 prokaryotic and 48 eukaryotic genomes were aligned and clustered on the basis of their sequence identity. This procedure led to the identification of classes of proteins that are peculiar to the different groups of organisms. Moreover the adopted similarity threshold produced clusters that are homogeneous on the structural point of view and that can be used for structural annotation of uncharacterized sequences. The second part of the work focuses on the characterization of thermostable proteins and on the development of tools able to predict the thermostability of a protein starting from its sequence. By means of Principal Component Analysis the codon composition of a non redundant database comprising 116 prokaryotic genomes has been analyzed and it has been showed that a cross genomic approach can allow the extraction of common determinants of thermostability at the genome level, leading to an overall accuracy in discriminating thermophilic coding sequences equal to 95%. This result outperform those obtained in previous studies. Moreover, we investigated the effect of multiple mutations on protein thermostability. This issue is of great importance in the field of protein engineering, since thermostable proteins are generally more suitable than their mesostable counterparts in technological applications. A Support Vector Machine based method has been trained to predict if a set of mutations can enhance the thermostability of a given protein sequence. The developed predictor achieves 88% accuracy.

Genomic characterization of the italian wolf (Canis lupus): the genes involved in black coat colour determination and application of microarray technique for snps detection.

This study provides a comprehensive genetic overview on the endangered Italian wolf population. In particular, it focuses on two research lines. On one hand, we focalised on melanism in wolf in order to isolate a mutation related with black coat colour in canids. With several reported black individuals (an exception at European level), the Italian wolf population constituted a challenging research field posing many unanswered questions. As found in North American wolf, we reported that melanism in the Italian population is caused by a different melanocortin pathway component, the K locus, in which a beta-defensin protein acts as an alternative ligand for the Mc1r. This research project was conducted in collaboration with Prof. Gregory Barsh, Department of Genetics and Paediatrics, Stanford University. On the other hand, we performed analysis on a high number of SNPs thanks to a customized Canine microarray useful to integrate or substitute the STR markers for genotyping individuals and detecting wolf-dog hybrids. Thanks to DNA microchip technology, we obtained an impressive amount of genetic data which provides a solid base for future functional genomic studies. This study was undertaken in collaboration with Prof. Robert K. Wayne, Department of Ecology and Evolutionary Biology, University of California, Los Angeles (UCLA).

Genomic and transcriptional analysis of allergen genes in apple (Malus x domestica)

Apple consumption is highly recomended for a healthy diet and is the most important fruit produced in temperate climate regions. Unfortunately, it is also one of the fruit that most ofthen provoks allergy in atopic patients and the only treatment available up to date for these apple allergic patients is the avoidance. Apple allergy is due to the presence of four major classes of allergens: Mal d 1 (PR-10/Bet v 1-like proteins), Mal d 2 (Thaumatine-like proteins), Mal d 3 (Lipid transfer protein) and Mal d 4 (profilin). In this work new advances in the characterization of apple allergen gene families have been reached using a multidisciplinary approach. First of all, a genomic approach was used for the characterization of the allergen gene families of Mal d 1 (task of Chapter 1), Mal d 2 and Mal d 4 (task of Chapter 5). In particular, in Chapter 1 the study of two large contiguos blocks of DNA sequences containing the Mal d 1 gene cluster on LG16 allowed to acquire many new findings on number and orientation of genes in the cluster, their physical distances, their regulatory sequences and the presence of other genes or pseudogenes in this genomic region. Three new members were discovered co-localizing with the other Mal d 1 genes of LG16 suggesting that the complexity of the genetic base of allergenicity will increase with new advances. Many retrotranspon elements were also retrieved in this cluster. Due to the developement of molecular markers on the two sequences, the anchoring of the physical and the genetic map of the region has been successfully achieved. Moreover, in Chapter 5 the existence of other loci for the Thaumatine-like protein family in apple (Mal d 2.03 on LG4 and Mal d 2.02 on LG17) respect the one reported up to now was demonstred for the first time. Also one new locus for profilins (Mal d 4.04) was mapped on LG2, close to the Mal d 4.02 locus, suggesting a cluster organization for this gene family, as is well reported for Mal d 1 family. Secondly, a methodological approach was used to set up an highly specific tool to discriminate and quantify the expression of each Mal d 1 allergen gene (task of Chapter 2). In aprticular, a set of 20 Mal d 1 gene specific primer pairs for the quantitative Real time PCR technique was validated and optimized. As a first application, this tool was used on leaves and fruit tissues of the cultivar Florina in order to identify the Mal d 1 allergen genes that are expressed in different tissues. The differential expression retrieved in this study revealed a tissue-specificity for some Mal d 1 genes: 10/20 Mal d 1 genes were expressed in fruits and, indeed, probably more involved in the allergic reactions; while 17/20 Mal d 1 genes were expressed in leaves challenged with the fungus Venturia inaequalis and therefore probably interesting in the study of the plant defense mechanism. In Chapter 3 the specific expression levels of the 10 Mal d 1 isoallergen genes, found to be expressed in fruits, were studied for the first time in skin and flesh of apples of different genotypes. A complex gene expression profile was obtained due to the high gene-, tissue- and genotype-variability. Despite this, Mal d 1.06A and Mal d 1.07 expression patterns resulted particularly associated with the degree of allergenicity of the different cultivars. They were not the most expressed Mal d 1 genes in apple but here it was hypotized a relevant importance in the determination of allergenicity for both qualitative and quantitative aspects of the Mal d 1 gene expression levels. In Chapter 4 a clear modulation for all the 17 PR-10 genes tested in young leaves of Florina after challenging with the fungus V. inaequalis have been reported but with a peculiar expression profile for each gene. Interestingly, all the Mal d 1 genes resulted up-regulated except Mal d 1.10 that was down-regulated after the challenging with the fungus. The differences in direction, timing and magnitude of induction seem to confirm the hypothesis of a subfunctionalization inside the gene family despite an high sequencce and structure similarity. Moreover, a modulation of PR-10 genes was showed both in compatible (Gala-V. inaequalis) and incompatible (Florina-V. inaequalis) interactions contribute to validate the hypothesis of an indirect role for at least some of these proteins in the induced defense responses. Finally, a certain modulation of PR-10 transcripts retrieved also in leaves treated with water confirm their abilty to respond also to abiotic stress. To conclude, the genomic approach used here allowed to create a comprehensive inventory of all the genes of allergen families, especially in the case of extended gene families like Mal d 1. This knowledge can be considered a basal prerequisite for many further studies. On the other hand, the specific transcriptional approach make it possible to evaluate the Mal d 1 genes behavior on different samples and conditions and therefore, to speculate on their involvement on apple allergenicity process. Considering the double nature of Mal d 1 proteins, as apple allergens and as PR-10 proteins, the gene expression analysis upon the attack of the fungus created the base for unravel the Mal d 1 biological functions. In particular, the knowledge acquired in this work about the PR-10 genes putatively more involved in the specific Malus-V. inaequalis interaction will be helpful, in the future, to drive the apple breeding for hypo-allergenicity genotype without compromise the mechanism of response of the plants to stress conditions. For the future, the survey of the differences in allergenicity among cultivars has to be be thorough including other genotypes and allergic patients in the tests. After this, the allelic diversity analysis with the high and low allergenic cultivars on all the allergen genes, in particular on the ones with transcription levels correlated to allergencity, will provide the genetic background of the low ones. This step from genes to alleles will allow the develop of molecular markers for them that might be used to effectively addressed the apple breeding for hypo-allergenicity. Another important step forward for the study of apple allergens will be the use of a specific proteomic approach since apple allergy is a multifactor-determined disease and only an interdisciplinary and integrated approach can be effective for its prevention and treatment.

Approcci molecolari e bioinformatici volti alla caratterizzazione di Vgt1, QTL coinvolto nella regolazione dell'epoca di fioritura in Zea Mays

The genetic control of flowering time has been addressed by many quantitative trait locus (QTL) studies. A survey of the results from 29 independent studies reporting information on 441 QTLs led to the production of a QTL consensus map, which enabled the identification of 59 chromosome regions distributed on all chromosomes and shown to be frequently involved in the genetic control of flowering time and related traits. One of the major QTLs for flowering time, the Vegetative to generative transition 1 (Vgt1) locus , corresponds to an upstream (70 kb) non-coding regulatory element of ZmRap2.7, a repressor of flowering. A transposon (MITE) insertion was identified as a major allelic difference within Vgt1. One of the hypotheses is that Vgt1 might function by modifying ZmRap2.7 chromatin through an epigenetic mechanism. Therefore, the methylation state at Vgt1 was investigated using an approach that combines digestion with McrBc, an endonuclease that acts upon methylated DNA, and quantitative PCR. The analyses were performed on genomic DNA from leaves of six different maize lines at four stages of development. The results showed a trend of reduction of methylation from the first to the last stage with the exception of a short genomic region flanking the MITE insertion, which showed a constant and very dense methylation throughout leaf development and for both alleles. Preliminary results from bisulfite sequencing of a small portion of Vgt1 revealed differential methylation of a single cytosine residue between the two alleles. ZmRap2.7 expression was assayed in the four developmental stages afore mentioned for the six genotypes, in order to establish a link between methylation at Vgt1 and ZmRap2.7 transcription. To assess the role of Vgt1 as a transcriptional enhancer, two reporter vectors for stable transformation of plants have been developed.

Wide-scale population genomics of Atlantic bluefin tuna (Thunnus thynnus) inferred by novel high-throughput technology

My PhD project was focused on Atlantic bluefin tuna, Thunnus thynnus, a fishery resource overexploited in the last decades. For a better management of stocks, it was necessary to improve scientific knowledge of this species and to develop novel tools to avoid collapse of this important commercial resource. To do this, we used new high throughput sequencing technologies, as Next Generation Sequencing (NGS), and markers linked to expressed genes, as SNPs (Single Nucleotide Polymorphisms). In this work we applied a combined approach: transcriptomic resources were used to build cDNA libreries from mRNA isolated by muscle, and genomic resources allowed to create a reference backbone for this species lacking of reference genome. All cDNA reads, obtained from mRNA, were mapped against this genome and, employing several bioinformatics tools and different restricted parameters, we achieved a set of contigs to detect SNPs. Once a final panel of 384 SNPs was developed, following the selection criteria, it was genotyped in 960 individuals of Atlantic bluefin tuna, including all size/age classes, from larvae to adults, collected from the entire range of the species. The analysis of obtained data was aimed to evaluate the genetic diversity and the population structure of Thunnus thynnus. We detect a low but significant signal of genetic differentiation among spawning samples, that can suggest the presence of three genetically separate reproduction areas. The adult samples resulted instead genetically undifferentiated between them and from the spawning populations, indicating a presence of panmictic population of adult bluefin tuna in the Mediterranean Sea, without different meta populations.

Functional and Genetic characterization of new genomic islands from an E. coli strain associated with neonatal meningitis.

Enterobacteriaceae genomes evolve through mutations, rearrangements and horizontal gene transfer (HGT). The latter evolutionary pathway works through the acquisition DNA (GEI) modules of foreign origin that enhances fitness of the host to a given environment. The genome of E. coli IHE3034, a strain isolated from a case of neonatal meningitis, has recently been sequenced and its subsequent sequence analysis has predicted 18 possible GEIs, of which: 8 have not been previously described, 5 fully meet the pathogenic island definition and at least 10 that seem to be of prophagic origin. In order to study the GEI distribution of our reference strain, we screened for the presence 18 GEIs a panel of 132 strains, representative of E. coli diversity. Also, using an inverse nested PCR approach we identified 9 GEI that can form an extrachromosomal circular intermediate (CI) and their respective attachment sites (att). Further, we set up a qPCR approach that allowed us to determine the excision rates of 5 genomic islands in different growth conditions. Four islands, specific for strains appertaining to the sequence type complex 95 (STC95), have been deleted in order to assess their function in a Dictyostelium discoideum grazing assays. Overall, the distribution data presented here indicate that 16 IHE3034 GEIs are more associated to the STC95 strains. Also the functional and genetic characterization has uncovered that GEI 13, 17 and 19 are involved in the resistance to phagocitation by Dictyostelium d thus suggesting a possible role in the adaptation of the pathogen during certain stages of infection.

Genomic and proteomic approaches in pig meat quality research field

Pig meat and carcass quality is a complex concept determined by environmental and genetic factors concurring to the phenotypic variation in qualitative characteristics of meat (fat content, tenderness, juiciness, flavor,etc). This thesis shows the results of different investigations to study and to analyze pig meat and carcass quality focusing mainly on genomic; moreover proteomic approach has been also used. The aim was to analyze data from association studies between genes considered as candidate and meat and carcass quality in different pig breeds. The approach was used to detect new SNP in genes functionally associated to the studied traits and to confirm as candidate other genes already known. Five polymorphisms (one new SNP in Calponin 1 gene and four additional polymorphism already known in other genes) were considered on chromosome 2 (SSC2). Calponin 1 (CNN1) was associated to the studied traits and furthermore the results reported confirmed the data already known for Lactate dehydrogenase A (LDHA), Low density lipoprotein receptor (LDLR), Myogenic differentiation 1 (MYOD1) e Ubiquitin-like 5 (UBL5), in Italian Large White pigs. Using an in silico search it was possible to detect on SSC2 a new SNP of Deoxyhypusine synthase (DHPS) gene partially overlapping with WD repeat domain 83 (WDR83) gene and significant for the meat pH variation in Italian Large White (ILW) pigs. Perilipin 1 (PLIN1) mapping on chromosome 7 and Perilipin 2 (PLIN2) mapping on chromosome 1 were studied and the results obtained in Duroc breed have shown significant associations with carcass traits. Moreover a study of protein composition of porcine LD muscle, indicated an effect of temperature treatment of carcass, on proteins of the sarcoplasmic fraction and in particular on PGM1 phosphorylation. Future studies on pig meat quality should be based on the integration of different experimental approaches (genomics, proteomics, transcriptomics, etc).

Metodi computazionali per l'annotazione di genomi e proteomi

Il progresso tecnologico nel campo della biologia molecolare, pone la comunità scientifica di fronte all’esigenza di dare un’interpretazione all’enormità di sequenze biologiche che a mano a mano vanno a costituire le banche dati, siano esse proteine o acidi nucleici. In questo contesto la bioinformatica gioca un ruolo di primaria importanza. Un nuovo livello di possibilità conoscitive è stato introdotto con le tecnologie di Next Generation Sequencing (NGS), per mezzo delle quali è possibile ottenere interi genomi o trascrittomi in poco tempo e con bassi costi. Tra le applicazioni del NGS più rilevanti ci sono senza dubbio quelle oncologiche che prevedono la caratterizzazione genomica di tessuti tumorali e lo sviluppo di nuovi approcci diagnostici e terapeutici per il trattamento del cancro. Con l’analisi NGS è possibile individuare il set completo di variazioni che esistono nel genoma tumorale come varianti a singolo nucleotide, riarrangiamenti cromosomici, inserzioni e delezioni. Va però sottolineato che le variazioni trovate nei geni vanno in ultima battuta osservate dal punto di vista degli effetti a livello delle proteine in quanto esse sono le responsabili più dirette dei fenotipi alterati riscontrabili nella cellula tumorale. L’expertise bioinformatica va quindi collocata sia a livello dell’analisi del dato prodotto per mezzo di NGS ma anche nelle fasi successive ove è necessario effettuare l’annotazione dei geni contenuti nel genoma sequenziato e delle relative strutture proteiche che da esso sono espresse, o, come nel caso dello studio mutazionale, la valutazione dell’effetto della variazione genomica. È in questo contesto che si colloca il lavoro presentato: da un lato lo sviluppo di metodologie computazionali per l’annotazione di sequenze proteiche e dall’altro la messa a punto di una pipeline di analisi di dati prodotti con tecnologie NGS in applicazioni oncologiche avente come scopo finale quello della individuazione e caratterizzazione delle mutazioni genetiche tumorali a livello proteico.

Next-Generation Sequencing-Based Mutations Scanning Strategy of the BCR-ABL Kinase Domain in Patients with PhiladelPhia-Chromosome Positive Leukemias Treated with Tyrosine Kinase Inhibitors

In chronic myeloid leukemia and Philadelphia-positive acute lymphoblastic leukemia patients resistant to tyrosine kinase inhibitors (TKIs), BCR-ABL kinase domain mutation status is an essential component of the therapeutic decision algorithm. The recent development of Ultra-Deep Sequencing approach (UDS) has opened the way to a more accurate characterization of the mutant clones surviving TKIs conjugating assay sensitivity and throughput. We decided to set-up and validated an UDS-based for BCR-ABL KD mutation screening in order to i) resolve qualitatively and quantitatively the complexity and the clonal structure of mutated populations surviving TKIs, ii) study the dynamic of expansion of mutated clones in relation to TKIs therapy, iii) assess whether UDS may allow more sensitive detection of emerging clones, harboring critical 2GTKIs-resistant mutations predicting for an impending relapse, earlier than SS. UDS was performed on a Roche GS Junior instrument, according to an amplicon sequencing design and protocol set up and validated in the framework of the IRON-II (Interlaboratory Robustness of Next-Generation Sequencing) International consortium.Samples from CML and Ph+ ALL patients who had developed resistance to one or multiple TKIs and collected at regular time-points during treatment were selected for this study. Our results indicate the technical feasibility, accuracy and robustness of our UDS-based BCR-ABL KD mutation screening approach. UDS was found to provide a more accurate picture of BCR-ABL KD mutation status, both in terms of presence/absence of mutations and in terms of clonal complexity and showed that BCR-ABL KD mutations detected by SS are only the “tip of iceberg”. In addition UDS may reliably pick 2GTKIs-resistant mutations earlier than SS in a significantly greater proportion of patients.The enhanced sensitivity as well as the possibility to identify low level mutations point the UDS-based approach as an ideal alternative to conventional sequencing for BCR-ABL KD mutation screening in TKIs-resistant Ph+ leukemia patients

Identification of novel genetic alterations in pediatric cytogenetically normal acute myeloid leukemia by next-generation sequencing

Pediatric acute myeloid leukemia (AML) is a molecularly heterogeneous disease that arises from genetic alterations in pathways that regulate self-renewal and myeloid differentiation. While the majority of patients carry recurrent chromosomal translocations, almost 20% of childhood AML do not show any recognizable cytogenetic alteration and are defined as cytogenetically normal (CN)-AML. CN-AML patients have always showed a great variability in response to therapy and overall outcome, underlining the presence of unknown genetic changes, not detectable by conventional analyses, but relevant for pathogenesis, and outcome of AML. The development of novel genome-wide techniques such as next-generation sequencing, have tremendously improved our ability to interrogate the cancer genome. Based on this background, the aim of this research study was to investigate the mutational landscape of pediatric CN-AML patients negative for all the currently known somatic mutations reported in AML through whole-transcriptome sequencing (RNA-seq). RNA-seq performed on diagnostic leukemic blasts from 19 pediatric CN-AML cases revealed a considerable incidence of cryptic chromosomal rearrangements, with the identification of 21 putative fusion genes. Several of the fusion genes that were identified in this study are recurrent and might have a prognostic and/or therapeutic relevance. A paradigm of that is the CBFA2T3-GLIS2 fusion, which has been demonstrated to be a common alteration in pediatric CN-AML, predicting poor outcome. Important findings have been also obtained in the identification of novel therapeutic targets. On one side, the identification of NUP98-JARID1A fusion suggests the use of disulfiram; on the other, here we describe alteration-activating tyrosine kinases, providing functional data supporting the use of tyrosine kinase inhibitors to specifically inhibit leukemia cells. This study provides new insights in the knowledge of genetic alterations underlying pediatric AML, defines novel prognostic markers and putative therapeutic targets, and prospectively ensures a correct risk stratification and risk-adapted therapy also for the “all-neg” AML subgroup.